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ABSTRACT 

We calculate the lowest-order non-linear contributions to the power spec- 
trum, two-point correlation function, and smoothed variance of the density field, 
for Gaussian initial conditions and scale-free initial power spectra, P(k) ~ k n . 
These results extend and in some cases correct previous work in the literature 
on cosmological perturbation theory. Comparing with the scaling behavior ob- 
served in N-body simulations, we find that the validity of non-linear perturbation 
theory depends strongly on the spectral index n. For n < —1, we find excellent 
agreement over scales where the variance u 2 {R) ^ 10; however, for n > —1, per- 
turbation theory predicts deviations from self-similar scaling (which increase with 
n) not seen in numerical simulations. This anomalous scaling suggests that the 
principal assumption underlying cosmological perturbation theory, that large- 
scale fields can be described perturbatively even when fluctuations are highly 
non-linear on small scales, breaks down beyond leading order for spectral indices 
n > — 1. For n < — 1, the power spectrum, variance, and correlation function in 
the scaling regime can be calculated using dimensional regularization. 
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1. Introduction 

Several independent arguments suggest that the growth of cosmological density per- 
turbations on large scales can be described by perturbation theory, even when the density 
and velocity fields are highly non-linear on small scales. For example, analytic and numer- 
ical work showed that linear perturbation theory describes the evolution of the large-scale 
density power spectrum P(k), provided the initial spectrum falls off less steeply than fc 4 
for small k QZerdovich 19651 , Peebles 1974] , Peebles fc Groth 197q , |Peebles 1980]) . In ad- 



dition, for Gaussian initial conditions, leading-order non-linear perturbative calculations of 
higher order moments of the density field, e.g., the skewness and kurtosis, agree well with 
N-body simulations in the weakly non-linear regime, where the variance of the smoothed 
density field o~ 2 (R) = (5 2 (R)) < 0.5 — 1 (Juszkiewicz, Bouchet & Colombi 1993, |Bernardeau 



|1994| , |Gaztahaga fc Baugh 1995| , Baugh, Gaztanaga & Efstathiou 1995). Thus, leading-order 
perturbation theory has been shown to work surprisingly well in the cases where comparisons 
with numerical simulations have been made. 

The success of leading-order cosmological perturbation theory raises questions: since 
the perturbation series is most likely asymptotic, what happens beyond leading order — does 
the agreement with simulations improve or deteriorate? More generally, what sets the limits 
of perturbation theory, beyond which it breaks down? These questions have become more 
urgent since it has been shown that leading-order perturbation theory appears to provide an 
adequate description even on scales where next-to-leading order and higher order perturbative 
contributions would be expected to become important. 

To address these questions, one must calculate loop corrections, i.e., corrections beyond 
leading order, in non-linear cosmological perturbation theory (NLCPT). This is the second 
paper of a series devoted to this topic. In the first paper ( [Scoccimarro fc Frieman 1996| ). 



we applied diagrammatic techniques to calculate loop corrections to 1-point cumulants of 
unsmoothed fields, such as the variance and skewness, for Gaussian initial conditions. Here, 
we calculate the 1-loop (first non-linear) corrections to the power spectrum, the volume- 
averaged two-point correlation function, and the variance of the smoothed density field for 
scale-free initial spectra, P{k) ~ k n . While the linear power spectrum for the Universe is 
not scale-free (on both observational and theoretical grounds), scale- free spectra are useful 
approximations over limited ranges of wavenumber k; they also have the advantage of yielding 
analytic closed form results. In a forthcoming paper, we will present 1-loop corrections to the 
bispectrum (the three-point cumulant in Fourier space) and the skewness of the smoothed 
density field for scale-free and cold dark matter spectra. 

One-loop corrections to the two-point correlation function and power spectrum have 
been previously studied in the literature flJuszkiewicz 1 981*1 , |Vishniac 1983| , Juszkiewicz, 
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Sonoda & Barrow 1984, Coles 1990| , Suto fc Sasaki 1991 , Makino, Sasaki & Suto 1992, 
Jain fc Bertschinger 1994] , [Baugh fc Efstathiou 1994] ). Multi-loop corrections to the power 
spectrum were considered by |Fry 1994] , including the full contributions up to 2 loops and 
the most important terms at large k in 3- and 4-loop order. Some of our results overlap in 
particular with the analytic results for the 1-loop power spectrum reported by Suto & Sasaki 
(1991) and in Makino et al. (1992). However, since we found that some of their expressions 
contain errors, we present complete corrected expressions for the 1-loop power spectrum. 
One-loop corrections to the variance were studied numerically for Gaussian smoothing by 
Lokas et al. (1995). We correct some numerical errors in their results, extend them to 
include top-hat smoothing and to the average two-point correlation function, and provide 
analytic derivation of some of the numerical results. 

The limiting behavior of the 1-loop corrections on large scales leads us to reconsider 
the issue of self-similarity in perturbation theory ( Pavis fc Peebles 1977 , Peebles 1980Q . 
It is commonly accepted that the evolution of density perturbations from scale-free initial 
conditions in an Einstein-de Sitter universe is statistically self-similar. N-body simulations 
have generally found self-similar scaling for —3 < n < 1, although the results for n < — 1 have 
been somewhat ambiguous due to problems of dynamic range in the simulations ( [Efstathiou 
et al. 1988] , [Bertschinger fc Gelb 199l| , [Ryden fc Gramann 199 lj , |Gramann 1992] , |Golombi~et 



al. 1995 , Jain 1995 , Padmanabhan et al. 1995|) . The issue of self-similarity in NLCPT has 
recently been investigated by Jain & Bertschinger (1995), who present arguments that the 
perturbative evolution is also self-similar. While we do not disagree with their calculations, 
we start from a different premise regarding cutoffs in the initial spectrum, which leads us 
to conclude that loop corrections break self-similar scaling if the spectral index n > — 1. In 
the regime where we do find scaling, n < — 1, we use dimensional regularization to calculate 
analytically the asymptotic behavior of the 1-loop power spectrum, variance, and average 
correlation function; these results agree quite well with the 'universal scaling' extracted from 
numerical simulations ( [Hamilton et al. 1991| , [Peacock fc Dodds 1994| , Jain, Mo & White 
1995). 

The paper is organized as follows. In Sec. § we discuss non-linear perturbation theory 
and review the diagrammatic approach developed in (|Fry 1984| , |Goroff et al. 1986| , |5cocci-| 
tnarro fc Frieman 1996|) to the calculation of statistical quantities. Section |3] presents the 



results of 1-loop calculations for the power spectrum, the smoothed variance, and average 
two-point correlation function. Self-similarity in perturbation theory is the subject of Sec. [|. 
We discuss the conditions under which 1-loop corrections exhibit self-similarity and compare 
our results with the universal scaling hypothesis, based on numerical simulations. Section |5] 
contains our conclusions. 
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2. Dynamics and Statistics 

2.1. The Equations of Motion 

Assuming the universe is dominated by pressureless dust (e.g., cold dark matter), in 
the single-stream approximation (prior to orbit crossing) one can adopt a fluid description 
of the cosmological N-body problem. In this limit, the relevant equations of motion cor- 
respond to conservation of mass and momentum and the Poisson equation (e.g., Peebles 
|1980| , [Scoccimarro fc Frieman 1996| ), 



dr 

<9v(x, r" 



(M(X " r) + V • {[1 + 5(x, r)]v(x, r)} = 0, (2-la) 



Or 



+ H{t) v(x, r) + [v(x, r) ■ V]v(x, r) = - V$(x, r), (2- lb) 

V 2 $(x, r) = ^ft 2 (r)<5(x, r) (2-lc) 



Here the density contrast 5(x, r) = p(x, r)/p — 1, with p(r) the mean density of matter, v = 
dx/dr represents the velocity field fluctuations about the Hubble flow, 7i = dlna/dr is the 
conformal expansion rate, a(r) is the cosmic scale factor, x denote comoving coordinates, r = 
/ dt/ais the conformal time, $ is the gravitational potential due to the density fluctuations, 
and the density parameter Q = p/ p c = 8irGpa 2 /3TL 2 . Note that we have implicitly assumed 
the Newtonian approximation to General Relativity, valid on scales less than the Hubble 
length aH~ x . We take the velocity field to be irrotational, so it can be completely described 
by its divergence 9 = V • v. Equations ( |2-1J) hold in an arbitrary homogeneous and isotropic 
background Universe which evolves according to the Friedmann equations; henceforth, for 
simplicity we assume an Einstein-de Sitter background, Q = 1 (with vanishing cosmological 
constant), for which a oc r 2 and 3Q7i 2 /2 = Q/t 2 . 

We take the divergence of Equation (|2-lfj| ) and Fourier transform the resulting equations 
of motion according to the convention 



for the Fourier transform of any field A(x, r), where, here and throughout, k is a comoving 
wavenumber. This yields 
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dS{k,r) 
dr 



0(k, r) = - / d 3 h / d 3 k 2 S D {k - k x - k 2 )a(k, ki)0(ki, r)5(k 2 , r), (2-3a) 



<9fl(k,r) 
dr 



H(r)6(k,T) + -n 2 (r)~5(k J r) = 

d 3 h f d 3 k 2 5 D (k - kx - k 2 )/3(k, k 1; k 2 )0(k 1; r)0(k 2 , r), (2-3b) 



(5d denotes the three-dimensional Dirac delta distribution), where the functions 



/, , x k • ki . . _ A; 2 (ki ■ k 2 ) 

a(k, ki) = — ^— , p(k, ki, k 2 J 



1.21.2 



;2-4) 



encode the non-linearity of the evolution (mode coupling) and come from the non-linear 
terms in the continuity equation ( |2-la| ) and the Euler equation ( [2-1 b| ) respectively. 



2.2. Perturbation Theory Solutions 

Equations ( j2-3| ) are very difficult to solve in general, being coupled integro-differential 
equations with no small parameter. The perturbative approach to the problem is to consider 
perturbations about the linear solution, effectively treating the variance of the linear fluctua- 
tions as a small parameter. In this case, Eqs. ( |2-3| ) can be formally solved via a perturbative 
expansion, 

oo oo 

5{k,r) = £ a n {r)8 n {k), 0(k,r) = H(r) £ a n {r)9 n (k), (2-5) 

n=l n=l 

where only the fastest growing mode is taken into account. At small a, the series are 
dominated by their first terms, and since 9i(k) = —5i(k) from the continuity equation, 5i(k) 
completely characterizes the linear fluctuations. The equations of motion (|2-3|) determine 
5 n (k) and 9 n (k) in terms of the linear fluctuations, 

5 n (k) = J d\ ... J d 3 q n 5 D (k - qi - ... - q n )Fi s \ qi , q n )5 1 (q 1 ) . . . ^(qj, (2-6a) 
6 n (k) = -J d 3 Ql ... J d 3 q n 5 D (k - qi - ... - q n )GW( qil . . . , q n )5i(qi) . . . 6 1 (q n ), (2-6b) 
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where F^ s ' and G$ are symmetric homogeneous functions of the wave vectors { qi, . . . , q n } 
with degree zero. They are constructed from the fundamental mode coupling functions 
ot(k, ki) and (3(k, ki, k2) according to the recursion relations (n > 2, see Goroff et al. (1986), 
Jain & Bertschinger (1994), for a derivation): 



F n (qi, . . . , q n ) = ^ 
in 



i (2n + 3)(n- 1) 
+ 2/3(k,k 1 ,k 2 )G 



(2n + l)a(k, ki)F n _ m (q m+ i, . . . , q r , 



(2-7a) 



C n (qi, . . . , q„ 



y-y C m (qi, . . . , q„ 



(2n + 3)(n-l) 
+ 2n/3(k,k 1 ,k 2 )G 



3a(k, k!)F n _ m (q m+ i, . . . , q„ 



(2-7b) 



(where ki = q x + . . . + q m , k 2 = q m+i + . . . + q„, k = k : + k 2 , and F 1 = G 1 = 1) and the 
symmetrization procedure: 



F^iqi, ■ ■ ■ , q n ) = — XI ^n(q7r(l)) • • • j q.7r(n))) 



^(qi, • • • > qn) = — ^(q^i), • • • , q^H), 

72*. — 



(2-8a) 
(2-8b) 



where the sum is taken over all the permutations it of the set {1, . . . ,n}. For example, for 
n = 2we have: 



if } (qi,q 2 ) 



<^ s) (qi,q 2 ) 



5 1 q x ■ q 2 /ft ft\ 2 (qx ■ q ; 

7 2 ft ft \q 2 ft/ 7 q\q\ 

3 1 qi ■ q 2 /ft. ft\ 4 (qi ■ q ; 

7 2 ft ft \q 2 ft/ 7 q\q\ 



(2-9a) 
(2-9b) 



1988 



The perturbation theory kernels have the following properties flGoroff ct al. 1986 , Wise 
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1. As k = q x + . . . + q„ goes to zero, but the individual do not, oc k 2 . This is a 
consequence of momentum conservation in center of mass coordinates. 

2. As some of the arguments of or G$ get large but the total sum k = qi + . . . + q„ 
stays fixed, the kernels vanish in inverse square law. That is, for p ^> q iy we have: 

FW( qi , . . . , q„_ 2 , p, -p) « GW( qi , . . . , q n „ 2 , p, -p) oc fc 2 /p 2 . (2-10) 

3. If one of the arguments qj of F^ or Gffl goes to zero, there is an infrared divergence of 
the form cfa/q 2 . This comes from the infrared behavior of the mode coupling functions 
a(k, ki) and (3(k, ki, k 2 ). There are no infrared divergences as partial sums of several 
wavevectors go to zero. 



2.3. Diagrammatic Expansion of Statistical Quantities 

In this work we focus on the non-linear evolution of two-point cumulants of the density 
field, such as the power spectrum and the volume-average two-point correlation function, 
and their 1-point counterpart, the variance. These are defined respectively by: 



;*(k, r)5(k', t)) c = 5 D (k + k')P(k, r), (2-11) 
£(R,t) = J £(x,T)W(x)d 3 x = J P(k,r)W(kR)d 3 k, (2-12) 

a 2 (R,r) = j P(k,r) W 2 (kR) d 3 k. (2-13) 

Here £(x, r) = / P(k, r) exp(ik • x) d 3 k is the two-point correlation function, and W(x) is a 
window function, with Fourier transform W(kR) which we take to be either a top-hat (TH) 
or a Gaussian (G), 

3 r i 
Wth(u) = — sin(u) — -ucos(-u) , (2-14) 

U A L J 

W G (u) = exp(-M 2 / 2 )- (2-15) 

We are interested in calculating the non-linear evolution of these statistical quantities 
from Gaussian initial conditions in the weakly non-linear regime, <r(R) < 1. A systematic 
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framework for calculating correlations of cosmological fields in perturbation theory has been 
formulated using diagrammatic techniques ( porofl et al. 1986| , |Wise 1988| , |Scoccimarro fc 



Frieman 1996| ). In this approach, contributions to p-point cumulants of the density field 
come from connected diagrams with p external (solid) lines and r = p — l,p, . . . internal 
(dashed) lines. The perturbation expansion leads to a collection of diagrams at each order, 
the leading order being tree-diagrams, the next to leading order 1-loop diagrams and so on. In 
each diagram, external lines represent the spectral components of the fields we are interested 
in (e.g., <5(k, r)). Each internal line is labeled by a wave-vector that is integrated over, and 
represents a linear power spectrum Pi±(q, r). Vertices of order n (i.e., where n internal lines 
join) represent a n th order perturbative solution 5 n , and momentum conservation is imposed 
at each vertex. Figure [I] shows the factors associated with vertices and internal lines. To 
find the contribution of order 2r to the p-point spectrum of the density field proceed as 
follows ( |Scoccimarro fc Frieman 1996| ): 



• Draw all distinct connected diagrams containing p vertices (with external lines labeled 
by ki . . . k p ) joined by r internal lines. Two diagrams are distinct if they cannot be 
deformed into each other by moving the vertices and lines without cutting any internal 
lines (sliding lines over other lines is allowed in the rearrangement process). For each 
of these diagrams: 

1. Assign a factor of <5o(ki — qi — ... — qn)^^ (qi> • • • , <in) to each vertex of order n 
and external momentum kj (i = 1, . . . ,p). For the arguments of F*f\ we use the 
convention of assigning a positive sign to wave-vectors outgoing from the vertex. 

2. Assign a factor of Pu(qj, r) to each internal line labeled by q,,. 

3. Integrate over all q., (j = 1, . . . , r). 

4. Multiply by the symmetry factor of the graph, which is the number of permuta- 
tions of linear fluctuations (Si's ) that leaves the graph invariant. 

5. Sum over distinct labelings of external lines, thus generating p\/{n\\ . . . n 2r _ p+ i!) 
diagrams (where rij denotes the number of vertices of order i). 

• Add up the resulting expressions for all these diagrams. 



To calculate 1-point cumulants, we have to integrate further over the kj (i = 1, . . .p). 
These integrations are trivial because of the presence of delta functions given by rule 1, so 
we are left only with integrations over the q/s. Also, since the p\/{ni\ . . .ri2r-p+i0 diagrams 
generated by rule 5 become equal contributions when the integration over external lines 
is performed, to find the contribution of order 2r to the p th -order 1-point cumulant of the 
density field we replace rule 5 by the following: 



• 5a. Integrate over kj (i = 1, . . . p) and multiply by the multinomial weight p\ / {n\\ . . . n^r-i 

Furthermore, when calculating cumulants of the smoothed density field, we add a window 
function at each vertex, giving instead of rule 1 (see Fig. |]): 

• la. Assign a factor of 5o(kj — qi — ... — q n )F^ s )(qi, . . . , q n )W(kiR) to each vertex of 
order n and external momentum kj (i = 1, . . . ,p). 

According to the above diagrammatic rules, we can write the loop expansion for the 
power spectrum up to one-loop corrections as 

P(k, r) = P i0 \k, t) + P (1) {k, r) + . . . , (2-16) 

where the superscript (n) denotes an n-loop contribution, the tree-level (0-loop) contribution 
is just the linear spectrum, 

P ( -°\k,r) = P 11 (k,r), (2-17) 

with a 2 (r)(5i(k)5i(k')) c = 5u(k + k')P 11 (/c, r), and the 1-loop contribution consists of two 
terms, 

P®(k, r) = P 22 (k, r) + P 13 (k, r), (2-18) 

with (see Fig. ||): 



•% . = S D (k - qx - . . . - q„)Fi s )( qi , . . . , q n )W(kR) 



q« v 



.2 = P u ( q , T ) 



Fig. 1. — Diagrammatic rules for vertices and internal lines for density field fluctuations 
smoothed with window W at radius R. 
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P 22 (k,r) = 2|[F 2 {s) (k-q,q)] 2 P 11 (|k-q|,r)P 11 (g,r)rf 3 g, (2-19a) 
P 13 (k,r) = 6j Ft\k,q,-q)P u (k,r)P 11 (q,r)d 3 q. (2-19b) 



Here Pij denotes the amplitude given by the above rules for a connected diagram representing 
the contribution from (SiSj) c to the power spectrum. We have assumed Gaussian initial 
conditions, for which P^ vanishes if % + j is odd. 

For the smoothed variance and average two-point correlation function we write 



a 2 (R) = a 2 e (R)(l + s^ a 2 e (R) + . . . 



(2-20a) 
(2-20b) 



where u 2 (R) and £e(R) denote the variance and average two-point correlation function in 
linear theory (given by ( |2-13 ) and ( [2-12 ) with P = Pn); the dimensionless 1-loop amplitudes 



are 



S M(R) 
x {1) (R) 



aj(R) 
1 

Wr) 



P (1) (A;,r) W 2 {kR) d 3 k, 
pM(k,r) W{kR) d 3 k. 



(2-21a) 
(2-21b) 



q 



1 1 k 



(Pi 



22) 



(P 



13 



Fig. 2. — Diagrams for the power spectrum up to one loop. See Eqs. ( |2-19 ) for diagram 
amplitudes. 
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A convenient property of these dimensionless amplitudes is the following: if one defines 
the correlation length Rq in linear theory as the scale where the smoothed linear variance 
is unity, a}{R^) = 1, then s^(Rq) = <j 2 (Rq^) — 1 is just the non-linear correction to the 
variance at this scale. 



3. Results 

3.1. One-Loop Power Spectrum 

We consider a linear (tree-level) power spectrum Pii(k,r) given by a truncated power- 
law, 

Pll(A; ' r) = \0 otherwise, ^ 

where A is a normalization constant, and the infrared and ultraviolet cutoffs e and k c are 
imposed in order to regularize the required radial integrations ( |Scoccimarro fc Frieman 1996| ) . 
In a cosmological N-body simulation, they would correspond roughly to the inverse comoving 
box size and lattice spacing (or interparticle separation) respectively. In the absence of the 
cutoffs, the spectrum (|3-1|) would be scale-free. We write the one-loop power spectrum 
contributions as 

Pij(k, r) = 7iA 2 a\r) k 2 c n+3 Pij (n; x, A), (3-2) 

where the dimensionless 1-loop spectrum = P22 + P13, and we have introduced dimen- 
sionless wavenumber variables, 

x ee k/k c , A ee k c /e. (3-3) 

From Eq. ( |2-19a|) we obtain (defining t = q/k c ): 



p 22 (n; x, A) ee — / dtt n dX (x 2 + t 2 - 2xt\) n/2 ~ 2 (3t - 7x\ - 10tA 2 ) 2 , (3-4) 

49 Ji/A Jx min (x,t) 
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when x < 2; otherwise P22{n;x, A) = 0. The constraint on the angular integration variable 
A = (k • q)/ (kq) comes from the cutoff dependence of -Pn(|k — q|, r), which gives: 



A m m(x, t) = Max < -1, — \ , 



(3-5a) 



f x 2 + t 2 _ A -2 ' 
X ma x{x,t, A) = Min U, - 



2xt 



(3-5b) 



Care must be taken when dealing with the limits imposed by Eqs. Q3-5a|) and (|3-5b| ), es- 
pecially in the cases n = —1, —2. On the other hand, integration over angular variables in 
Eq. ( |2-19b|) is straightforward, and we obtain: 



p 13 (n;x,A) = x r ' 



x In 



dt t n+2 

1/A 
X + t 



(Qx 6 - 79xH 2 + 50x 2 t 4 - 21t 6 ) (t 



x 



2\3 



(7t 2 + 2x 2 



63x 2 t 4 



A2xH 5 



\x — 1\ 



(3-6) 



for A 1 < x < 1; otherwise Pi3(n; x, A) = 0. 

In the following subsections we give results for pi% and P22 in the limit A = k c / e — ► 00 
up to terms of order A , for spectral indices n = 1, 0, —1, —2. We also give their asymptotic 
behavior at large scales. From the properties of the perturbation theory kernels given in 
Sec. |2]2| and Eqs. (|2-19|), one would naively expect that as x — ► 0: 



P 22(n; x, A) oc x 4 ~ (k/k c )\ p 13 (n; x, A) oc x n+2 ~ (k/k c ) n+2 . (3-7) 

Although this is certainly correct for p 2 2 when x A -1 (i.e., k e), the expressions below 
correspond to P22 and P13 after the limit A — > 00 has been taken and therefore exhibit a 
different kind of asymptotic behavior (corresponding to <— e < k k c ). This deviation 
from the expected scaling becomes more pronounced as n decreases, since infrared effects 
(e — > 0) become more important with the increase of large-scale power. As we will see, 
Eq. ( |3-7[ ) is obeyed by P22 only when n > 1 and by piz when n > 0. 

Our results below are equivalent to those of Makino et al. (1992) in the case of pis, 
but they differ for p 2 2 hi two respects. First, their expressions for p 2 2 are only valid for 
x < 1; they did not seem to consider that P22 is non- vanishing up to x = 2, and the region 
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1 < x < 2 requires a separate integration when n is odd. In addition, for n even, their series 
expansions for the dilogarithms do not converge in this region. The second difference comes 
from the fact that for n = —1,-2, p 22 develops a divergence at x — 1, and the expressions 
for P22 must be modified in a region of radius A -1 about x = 1. This singularity is integrable, 
giving a finite contribution to the 1-loop correction to the variance. Finally, some of their 
resulting expressions for p^ contain typographic errors. We have checked our expressions by 
analytically integrating them and verifying that they correctly reproduce the unsmoothed 
1-loop coefficients given by Scoccimarro & Frieman (1996). For comparison, we also 
plot (but do not give explicit expressions for) the 1-loop power spectrum in the Zel'dovich 
approximation. 



3.1.1. n = l 

For n=l, the dimensionless 1-loop power spectrum contributions are 



, , 1 52a; 181a; 3 2x 5 , /l + x\ (5 - 19x 2 + 25x 4 - 5a; 6 + 10a; 8 
p 13 {l;x) = -77— + — — — + ln' 



21a; 315 315 21 VI -x) 210x 2 
8x 5 , /l — x 2 \ . 

"Ws h <— )• (3 - 8) 



, , 18x 4 13x 5 20x 6 x 7 . . 

p 22 (l ]X ) = + — , ifx<l, (3-9) 

v V ; 49 98 1029 49' ~ ' V ; 



, s 320 I62 2x 3 18a; 4 5x 5 20x 6 x 7 , 

p 22 (l-x) = + + + + — , if 1< x < 2, 3-10 

y K 1 J 1029x 49 3 49 42 1029 49 _ _ , v ; 

and the 1-loop reduced power spectrum is then (see Fig. 131) : 



m , , /, ^ 1 52x 181x 3 18x 4 67x 5 20x 6 x 7 

p w (n = l:x = k k c < 1) = 1 1 1 

F y 21x 315 315 49 294 1029 49 

1 + x\ (5 - 19a; 2 + 25x 4 - 5x 6 + 10a; 8 ) 
210a; 2 

(3-11) 




-14- 



with p^(x) = p 2 2(x) when 1 < x < 2. For i< 1 we get: 




(3-12) 
(3-13) 



J i i i i I i i i i J i i i i I i i i i L 

0.5 1 1.5 2 

x 

Fig. 3. — One-loop corrections to the power spectrum for n = 1. The solid line shows 
p (1) (l;ar) as a function of x = k/k c . Dashed lines correspond to the contribution of individual 
diagrams, p 2 2(l;^) (long-dashed) and pi 3 (l;x) (short-dashed), whereas the dot-dashed line 
shows pW(l;x) in the Zel'dovich approximation. 



3.1.2. n = 

For n = 0, the 1-loop power spectrum contributions are 



, , 1 157 2Q9x 2 n 2 x 3 x 4 , (\ 

P1 ^ 0;X) = -T8^ + 756 -^52--^r-21 +ln ll 



x 



x^— ^(-14 + 43a; 2 - 47x 4 + 12a; 6 ) + ^[Li 2 (a;) - Li 2 (-a;)], (3-14) 
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p 22 (0;x<2) 



25 25x 

98 + 196 ~ 147~ ~~ 392 



10x 2 5x 3 29vr 2 x 3 3x 4 65x 5 
+ 



294 + 49 



1176 196 



+ In | x — 1 1 



X 



(25 - 15x 2 - 16x 4 + Qx 6 + 29x 4 lnx) 29x 3 rl 



98a; 



98 



ln 2 (a;) + Re[Li 2 (x)] 



+Li 2 



x — 1 



X 



(3-15) 



where Li 2 denotes the dilogarithm, defined by ( [Lewin 1981| ) 



Li 2 (a;) 



ln(l - z) 



dz 

o z 



(3-16) 



which has the series expansion for small argument [x < 1): 



m=l 



2 ' 



(3-17) 



The resulting 1-loop power spectrum is (Fig. |j) 



1 2449 25x 2003x 2 5x 3 17vr 2 x 3 2x 4 65x 5 x 7 



18x 2 5292 196 1764 392 196 147 1176 196 

+ in(l + X )S± M-U + 43x 2 - 47x 4 + 12x 6 ) - — l n 2 (x) 

v ; 504x 3 V ; 196 V ; 

(-98 + 1299x 2 - 1170x 4 - 163x 6 + 132x 8 + 1044x 6 In x) , , 

3528x d 

_!^ L i 2 (x) - — Li 2 (-z) - — Li 2 ( , (3-18) 
147 V ; 42 2V ; 98 V x J ' V ; 



and p^'(x) = p 22 (x) when 1 < x < 2. 
For x< 1 we have: 



/ \ 244 

Pi 3 (0;x) « P {1) (0;x)^-— x 2 + 0(x 3 ), (3-19) 

315 

29tt 2 

P22(0;x) « -— x 3 + 0(x 4 ). (3-20) 
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3.1.3. n= -1 

For n < — 1, the 1-loop contributions are infrared-divergent and therefore depend on 
the infrared cutoff through A = k c / e, 



4i, , 1 88 lb 2x 3 , fl + x 
p 13 (-l;x,A) = -—InA--— + ——- — - — + In 



3 15s 3 315s 189 63 \l-x, 

(21 - 95x 2 + 225x 4 + 15x 6 + 10x 8 ) 88x , / 1 - x 2 \ 



p 22 {-l;x,A) 



Ax 



lnA4 
3 147 



80x 2x 2 



x 



+ 



44x 4 2x r ° 
4 



441 49 441 



+ 



2x 



ln[(l-x)x 2 ], 



x < l-A- 1 , 
(3-22) 



p 22 (-l;x, A) 



— InA 

3 



160 x 3 2x 5 x 7 2i, , 

1 1 1 1 ln(x 

147x 3 49 441 3 v 

.2 , i n„,3 i c„4\ a ( X ~ 



4-105x 2 4- 12x 3 4- 6x 4 ) A 



882x 



(x - \\ 2 

v — ; (50 4- 100x 
147x 

(15 + 45x 4- 60x 2 4- 60x 3 4- 12x 4 
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Fig. 5. — Same as Fig. ||| but for n — — 1 and A = k c /e = 100. 



+4x & )A d , 1 - A" 1 < x < 1 + A 



-l 



(3-23) 



320 



80x 2x 2 x 3 44x 4 2x 5 x 7 



U7x 147 



The 1-loop correction is then 



3 + 3 441 



2x 



I 1 ln(x-l), 

49 441 3 V ; ' 



1 + A" 1 < x < 2, 
(3-24) 



V 



T;xl 



1 88 
+ 



643x 2x 2 19a; 3 44x 4 2x 5 x 7 

15x 3 ' 315a; + 1323 + ~3~ + "63" + HT + 49" + 441 

(21 - 95x 2 + 225a; 4 + 15a; 6 + 10x 8 ) , / 1 + x \ 298a; , , 
+- 7^7— > In ( ) + ln(l - x) 



630a; 4 

244x, . . 88x Wl , 

in x H m l + x) 

315 V ; 315 V ; 



. 1 — x 
x < 1 - A -1 , 



315 



(3-25) 



whereas when 1 < x < 2 we have 1; x) = P22(— 1; sc)- For a; C 1 we get: 



p 13 (-l;x, A) 



4x, . 128 176x, , , ^ 

In A H x ln(x) + 0(x 3 ) 

3 225 315 V ; V ; 



(3-26) 
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p 22 (-l;x,A) 
p«(-l;x) 



y In A + —x + — \n(x) + 0(x 4 ), 

12272 244a;, , . , nl 

: x + -^rln(a;) + 0(aT). 



11025 



315 



(3-27) 
(3-28) 



In this limit, the divergences cancel; note, however, that the A — > oo divergence in p 22 in 
the region 1 — A -1 < x < 1 + A -1 is uncancelled, although the 'size' of the divergent region 
shrinks to zero. 



II.4. n 



For n = —2, the 1-loop terms are 



Pl3 (-2;x,A) 



- 4 -A- 

3 12a; 4 



1 107 5vr 2 82 x 2 , 
H r H 1 h In 

252a; 2 28a; 63 42 VI 

5 



1 +x\ (x 2 - 1) 



x 



-7 + 31x z + 4a; 4 + 2x D + 



14a; 



[Li 2 ( 



X 

Li 2 (a;)], 



168a; 5 



(3-29) 



p 22 (-2;x,A) 



2 . 75 
-AH 

3 1568a; 



25tt 2 15x 29a; 3 



196a; 392 
-29a; 4 - 3a; 6 - 30a; 2 In a;) + 



392 
A 4 



3a; 5 
196 



x 

784 



A 5 



196a; 



-50 + 82a; 2 



+44a; 4 lna;), 1 



1568a; 
A- 1 <x < 1 + A^ 1 



(25 - 60a; 2 + 63a; 4 - 24a; 6 - 4a; 8 



(3-30) 



whereas when a; < 1 — A 1 or a; > 1 + A 1 we have: 



p 22 (-2;x,A) = 



4 1 

-e(i-x)A + - 

3 V ; 2(x-l) 



25n 2 



205 309a; 11a; 2 



53a; 3 



x 



75 



784 
75 



98a; 



98a; 

ln 2 (a;) + In \x - 1 
x — 1 



392 784 392 1568 
;-1001 + 60a; 2 - 11a; 4 + 300 In a;) 



3a; 5 
196 



196a; 
Re[Li 2 (x)] + Li 2 



392a; 



x 



(3-31) 



where 6{x) is the step function. This gives: 
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(1) 1 107 1 85vr 2 2747 309x 61x 2 53x 3 3x 5 

P ^~ ' X ' ~ 12x 4 + 252x 2 + 2(x - 1) + 196x + 3528 ~ 784 ~ 1176 + 1568 + 196 

x 7 75 , ,. . , , . (-1001 + 60x 2 -llx 4 + 300 lnx) 
+ — - — — ln 2 (x) + ln(l - x) v 



784 196x v ' y ' 392x 
+ In [\^\ ( f "^ (-7 + 31x 2 + 4x 4 + 2x 6 ) + ^-U 2 (-x) 



1 — xj 168x 5 14x 
55 Li 2 (x)-^Li 2 (^), x<l-A-\ (3-32) 



49x 98x \ x 

otherwise 2; x) = p 22 (— 2; x) when 1 + A -1 < x < 2. For x < 1 we have: 



p 13 {-2;x,A) « -!a + ^ + O(x ), (3-33) 
4 757r 2 

p 22 (-2;*,A) » A+ +O(x ), (3-34) 
3 196x 

55-7T 2 

p«(-2;x) « +0(a; ). (3-35) 
98x 



3.2. One-Loop Smoothed Variance and Average Two-point Correlation 

Function 

Using the results of the previous section for the power spectrum, we can calculate the 1- 
loop corrections to the variance and average two-point function, Eqs. (|2-13| , 2-12 ). We focus 
on the corresponding 1-loop amplitudes and x^ defined in Eqs. ( 2-20 , 2-21|) . From 
( ggg ), Q , and O we can write: 



s {1 \n;k c R) e 
x {1) (n; k c R) 



(k c R) 



2n+3 



2k c R 



4[A/ CT (n, k c R,A)]2Jo 

{k c R) 2n+3 r2k c R 



4[AI e (n, fc CJ R, A)] 2 Jo 



p (1 \n; u(kcR)^ 1 , A) PU 2 (w) w 2 dw, (3-36a) 
pW(n;u(A; .R)-\A) W(u) w 2 du, (3-36b) 



where AI a (n, k c R, A) = I CT (n, fc c i?) — I a (n, eR) and similarly for A/g, with 
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Un,z) = / u n+2 W 2 (u) du, (3-37a) 
Jo 

Iz(n,z) = j Z u n+2 W{u) du. (3-37b) 



o 



In Appendix [A] we present analytic results for I a and 1^ with top-hat and Gaussian smoothing. 

In Figs. [7| - we show results for and x^ for top-hat and Gaussian smoothing, 
based on numerical integration of Eq. (|3-36|) . Similar results for s^ 1 ' for the Gaussian window 
function were presented by Lokas et al. (1995); comparison with our results shows very good 
agreement for scales such that k c R > 2. However, in the limit k c R — > 0, they apparently find 

~ 0, instead of recovering the unsmoothed values reported in Scoccimarro & Frieman 
(1996) and in conflict with our results shown here. The source of their error may have been 



to replace 2k c R by k c R in the upper limit of integration in Eq. ( |3-36a| ); we found we could 
approximately reproduce their results with this change. 

From Eqs. ( |3-36|) , one can analytically calculate sW and in the limit k c R ^ 1 by 
asymptotic expansion in 1/ k c R: for large k c R, all we need is the behavior of p^ (n; x, A) for 
small x given in the previous section. For n = —2 we have: 

S W(-2; k c R » 1) « s W(-2; oo) = — r °°) , (3-38) 

V ' ; v ' ; 392 [J (J (-2,oo)] 2 ' V ; 



where we used Eqs. (|3-35|) and (|3-37a|) . An identical expression holds for x^ 1 ' upon replacing 



io-'s by I^s. Using the results in Appendix |A], we find for top-hat smoothing: 



1 375 

sS(-2;oo) = — ^ 0.877, (3-39) 
4l(-2;oo) = ^« 0.748. (3-40) 



For Gaussian smoothing we find 

s£\-2; oo) = sg>(-2; oo) = || « 0.882 (3-41) 

in good agreement with the numerical results of Lokas et al. (1995), who found Sq\—2; oo) ~ 
0.86. Comparing these analytic results to Fig. [7], we see excellent agreement except for x^y, 
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which oscillates around the value given by Eq. ( |3-40|) . The reason for these oscillations can 
be partially understood from Eq. ( |3-38| ), using the results in Appendix |A] for finite k c R. The 
behavior of x^e. is a l so due in part to the fact that smoothing with one top-hat window 
allows small-scale power in p^(n; x, A) to 'leak in' to the average correlation function, while 
( |3-38|) assumes that only the large-scale power contributes. For this reason, in the following 
we will give analytical results only for Gaussian smoothing. 

For n = — 1, using Eqs. ( |3-28|) and ( [3-3 7a|) we have: 



ft!/ i i n ^ 3068 J - 1 1 ' 00 ) 244 r°°Ti,2/ \ 3 1 7 61 J ff 1,00 

s (j (-l; k c R > 1) = F P — ^— H / jy 2 (u)w 3 lnw r P — 

V ' ; 11025 [I CT (-1, oo)] 2 315 Jo K J 315 [J ff (-1, oo)] 2 

x \n(k c R) (3-42) 
and similarly for Using Gaussian smoothing we find: 



4 1) (-1;A; CJ R»1) = + |i(l - 7e ) - ^ ln(fc cJ R) « 0.638 - 0.387 \n(k c R), (3-43) 

6136 61 122 
4 1} (-l;A; c i?>l) = __ + _(l- 7e + In 2)-— ln(A: CJ R) w 0.773 - 0.387 ln(A; c i2) 

11025 315 315 

(3-44) 



where 7e ~ 0.577216... is the Euler-Mascheroni constant. Eqs. (|3-43|) and (|3-44j) are plotted 



as the dot-dashed curves in Fig. |8|, which shows the excellent agreement with the results 
from the full numerical integration at large k c R. 

For n = 0, one must include the next to leading order term at small x in p^(0; x, A), 
since this gives rise to a constant term in the 1-loop coefficients. We have: 

P {1 \0;x) « -^* 2 + ^* 3 + 0(x 4 ), (3-45) 
which gives (see dot-dashed curves in Fig. ||): 



4! } (0;£; CJ R> 1) = — k c R « 1.710 - 0.655 A; CJ R, (3-46) 

; 147 105v^F c , v ; 

x[! } (0;A; CJ R> 1) = - — J- k c R « 1.710 - 0.464 A; c i2. (3-47) 

; 147 105 V 7T v ' 
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Similarly, for n = 1, we need the expansion: 

p u> (1; „ M _ 1», + ^ -(^ + ^ln I )/ + 0( A (3-48) 
which leads to (see Fig. |i~0l): 



4 1) (l;fc cJ R»l) = -^-1(11-^)1^(^)1^^-5^ 
GV ' ; 7350 105 V [ ; 35 V ; 392 315 V ; 

w -0.964 + 0.457 \n(k c R) + 0.610 k c R - 0.387 (A; CJ R) 2 , (3-49) 

xWtl-.kcR^l) = -l(ll-67 e + 61n2) + — \n(k c R) + — \[^k c R 

G v ' ; 7350 105 V ' ' 35 V ; 392 V 2 

61 

-(fc c i?) 2 ss -1.122 + 0.457 \n(k c R) + 0.432 k c R - 0.194 (fc c i?) 2 . 



315 

(3-50) 

Note that the 1-loop coefficients are scale-dependent for n > — 1; this is an indication of 
the breaking of self-similarity, as we discuss next. (A similar analysis to that above can be 
carried out for the 1-loop corrections in the Zel'dovich approximation, but details will not 
be given here.) 



4. Self-Similarity and Perturbation Theory 

4.1. Self-Similar Solutions 

Since there is no preferred scale in the dynamics of a self-gravitating pressureless perfect 
fluid in an Einstein-de Sitter universe, Eqns. (|2-1|) admit self-similar solutions (see |Peebles| 
1980|) . This means that the cosmological fields should scale with a self-similarity variable, 



given appropriate initial conditions: knowing the fields at a given time completely specifies 
their evolution. We can search for the appropriate self-similarity transformation by rewriting 
the fields as 



5(x,r) = x (y), (4-la) 
v(x,r) ee r"u(y), (4-lb) 
$(x,r) ee rVfr), ( 4 " lc ) 
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where y = x/t" is the similarity variable. Self-similarity is obeyed if we can rewrite Eqs. 
[I]) in terms of y only; this condition determines the indices /x, 7, and v. Note that the term 
1 + <5(x, r) in the continuity equation fl2-la| ) precludes a power of r prefactor in Eq. ( |4-la|) . 



Substituting the ansatz Eq. (|4-1|) into Eq. ( |2-la| ) and equating powers of r we get the 
condition 



H = v-1, (4-2) 



whereas the Poisson equation Eq. ( 2-lc ) gives: 



7 = 2(z/-l). (4-3) 

These two conditions in turn guarantee that the Euler equation obeys self-similar evolution. 
The index v can be determined in terms of the spectral index n of the linear power spectrum: 
Eq. ( |4-la|) implies that the two-point correlation function obeys the self-similar scaling 

C(x,r)=E(xr-n, (4-4) 
which upon Fourier transformation gives for the power spectrum 



P(k,T)=r 3u V{kr v ). (4-5) 

On the other hand, in the linear regime, for a power-law initial power spectrum with spectral 
index n in an Einstein-de Sitter universe, Pn(k, r) ~ r 4 fc n ; equating this with (|4-5|) fixes the 
remaining parameter in the similarity transformation, 

(A more general derivation of ( |4-(j| ) starts from the BBGKY hierarchy of equations for the 
phase space particle distribution functions, |Davis fc Peebles 1977| .) Using these results, the 



equations of motion (|2-1D can be rewritten in self-similar form: 



vy ■ Vx(y) - V • {[1 + x(y)]u(y)} = 0, (4-7a) 
uy ■ Vu(y) - (1 + u)u(y) - [u(y) ■ V]u(y) = V^(y), (4-7b) 



vV(y) = 6x(y), 



(4-7c) 



where V = d/dy. The solutions of these equations in the linear regime are x(y) ^ with 
( = — 2/u, 3/u, which correspond to the well-known linear growing and decaying modes, 
5i oc a,a~ 3 / 2 . The condition of self-similar evolution also fixes the spatial dependence of 
the fields, e.g., <5i(x) oc x~ 2 ^ u for the growing mode. However, we generally consider the 
density and velocity fields to be random fields, on which we impose initial conditions only 
statistically, e.g., by specifying the linear power spectrum for Gaussian initial conditions. A 
particular realization of the ensemble of initial conditions will not obey this spatial scaling, 
so the cosmological fields themselves will not be self-similar QJain fc Bertschinger 19951 ). 
Since NLCPT is built from linear solutions, the same conclusion applies to the expansions 
given by Eq. (|2-5|) , which are not self-similar (except when the spatial dependence is fixed as 
above). Nevertheless, in some cases the statistical quantities of interest, such as the power 
spectrum, variance, and two-point correlation function, may exhibit self-similar scaling even 
if the fields in a given realization do not. We now consider the conditions under which this 
happens in perturbation theory. 



4.2. Self-Similarity and Linear Perturbation Theory 

The introduction of fixed (time-independent) cutoff scales e and k c in the linear power 
spectrum ( |3-1| ) breaks self-similarity, because they do not scale with the self-similarity vari- 
able fca 2// ( n+3 ) (the Fourier space analog of y). The extent to which one can take the limits 
e — > and k c — > oo will determine whether the statistical properties of the density field obey 
self-similar scaling. In the absence of cutoffs, the only physical scale that can be defined 
from the power spectrum is the correlation length, Rq, defined by: 

a 2 (R ) = J P(k,T)W 2 {kR )d 3 k = 1. (4-8) 

In linear perturbation theory, for scale-free initial conditions, one finds Rq oc a 2 ^ n+3 \ which 
has the right time dependence to build the self-similarity variable. Consequently, in the 
absence of cutoffs, statistical quantities in linear theory evolve self-similarly with i?o, and we 
can write, e.g., (see Eqs. ( |4-4|) and (p5|)): 



Ce(r,r) = E £ (r/Ro), 



(4-9) 
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R^P u (k,T)=V n (kR ). (4-10) 

The presence of cutoffs in the initial spectrum changes this situation. From Eq. ( |3- 1\ ) 
we have 

rk c Ro 

R™ +3 = 4nAa 2 / K n+2 W 2 {n)dK = 4nAa 2 [I a (n, k c R ) - Un, eR )} 

JeRo 

= AnAa 2 AI a (n,k c R ,A), (4-11) 

and Rq will not scale as a 2 /( n+3 ) unless AI a (n, k c Ro, A) is time independent. Before examining 
under which conditions this is true, it is useful to introduce another measure £q of the 
correlation length, defined in terms of the volume-averaged correlation function, 

£(4) = J P{k, r)W(k£ )d 3 k = 1, (4-12) 

and therefore: 



= AnAa 2 f ° ° K n+2 W{K)dK = 4nAa 2 [I^n, kj ) - I^n, el Q )\ = 4nAa 2 AI^n, kj , A). 

(4-13) 



Note that Rq and £q both differ slightly from the conventional definition of the correlation 
length, tq, in terms of the two-point correlation function, £(ro) = 1. 

A necessary condition for self-similarity in the variable Rq (£q) is the convergence of the 
integrals AI a (AI^) for small and large n (where k denotes kR (k£ )); otherwise these quan- 
tities would be sensitive to the cutoffs and consequently time dependent. From Eqs. ( 4-ll| ) 
and ( |4-13| ), in order for the correlation length to scale self-similarly in linear theory, we must 
have at least the following conditions: 



1. For a Gaussian filter, since W(k) w 1 when k — > 0, convergence in the infrared (small 
k) requires n > —3. Convergence in the ultraviolet (large k) is achieved for any n. 

2. For a top-hat filter, convergence in the infrared requires n > —3. Since Wtk{k) ~ k~ 2 
as n — ► oo, convergence in the ultraviolet requires n < 1 in Eq. ( [1-1 1|) and n < — 1 in 
Eq. (EH- 
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Note the dependence of these conditions on the local averaging scheme (top-hat vs. Gaussian) 
and the type of statistics (variance vs. average correlation function). 

While these conditions guarantee convergence of AI a ^, they do not test the constancy of 
these quantities for different choices of parameters. In Figs. |ll]- [14], we show AI a ^(n, k c Ro, A) 
for A = k c /e = 10 2 , 10 3 , oo, as a function of k c R for spectral indices n = —2, —1, 0, 1. These 
finite values of A are comparable to the dynamic range currently achievable in N-body 
simulations. By definition, k c Ro measures the ratio of the correlation length to the small 
scale cutoff of the linear power spectrum. Therefore, in the region k c Ro « 1, the evolution of 
the correlation length is strongly affected by the cutoff k c , and the AI a ^ drop precipitously. 
By the time the correlation length has evolved to k c R > 3, however, AI a ^ reach their no- 
cutoff limits (in the cases where they converge). On the other hand, when k c R Q « A, i.e., 
R pa e" 1 , the correlation length reaches the large-scale cutoff, and AI a ^ again drops. This 
behavior is more dramatic as n decreases (and therefore strongest for n = —2) due to the 
increase in large-scale power. As expected from the discussion above, the definition of the 
correlation length using Gaussian smoothing displays self-similar evolution over a longer time 
interval than the other possibilities, with top-hat smoothing becoming notoriously worse as 
n increases. 

The behavior of statistical quantities in linear perturbation theory with respect to self- 
similarity is determined directly from these considerations. We can write the linear power 
spectrum as (see Eq. ( |4-10| ) ): 

V n (kR 0} k c R ,A) = ^ R f (4-14) 
AirAl [n, k c R , A) 

where AI depends on the local averaging scheme and statistical quantity used to define 
the correlation length, and kRo <C 1 in the linear regime. Unless otherwise noted, for power 
spectrum calculations we will use the choice in Eq. (|4-llj ), i.e., Rq, with Gaussian smoothing. 



For the variance and average correlation function we have 



^^^^M/iy 1 " 1 , (4 . 15) 



AI a (n } k c R 0} A) \R { 



»o 



-(n+3) 



Now the question is under what conditions are the statistics self-similar in the variables 
kRo (for the power spectrum) and R/Ro (variance and average correlation function). In the 



regime R 3> Rq (where linear theory is valid) and 



3 < A; c i? < k c R < A, (4-17) 



Figs. |TT| - |TJ show that AI is approximately constant. Provided the conditions on n stated 
above hold, in this regime kRo and R/Rq are self-similar variables: in linear perturbation 
theory, the power spectrum scales as 

V n (kR Q ,k c Ro,A) « ( ^ o) " , (4-18a) 
2nT (^) 



(see Eq. (TAjJ) ), and 



at(R/Ro, k c R , A) « &(i?/#o, Mfc, A) « ( — j . (4-18b) 



Note that, for top-hat smoothing, the requirement of self-similarity of the statistical quan- 
tities is in practice less restrictive than the requirement that Rq or £ scale self-similarly. 
For example, although defined from top-hat smoothing, requires n < — 1 to scale self- 
similarly, the average correlation function for n = 1 does scale self-similarly, i.e., it depends 
on R/Rq only, in the regime given by Eq. ( |4-17| ), 



UR/Ro,k c Ro,^) 



R_ 

Rq 



(4-19) 



(see Eq. ( |A8d| )) where we have averaged over oscillatory behavior. However, it does not obey 
the "expected" £ e (R/R ,k c R Q ,A) w (R/R )- 4 scaling. This trend appears to be seen in the 



results from high-resolution scale- free N-body simulations (Colombi et al. 1995). 



4.3. Self-Similarity and Non-Linear Perturbation Theory 

We have examined the conditions under which linear perturbation theory obeys self- 
similar scaling. We now turn to the question of whether loop corrections in NLCPT affect 
self-similar scaling. 

Using the results of the previous section and Eq. fl3-2| ) , we can write the power spectrum 
up to 1-loop corrections as: 
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v(kR , k cRo , a) = A /^°, ) " . \ + 1fi r A ( ^f o) rr a^ 12 p {i) ^ k i k - ( 4 - 2 °) 

47rA/ (T (n, fe c ito, A) lo7r|Ai CT (n, k c R , A)\ z 

where kRo < 1 in the weakly non-linear regime. Self-similarity is maintained at the 1-loop 
level if, in the scaling regime given by 



kR < 1 < k c R < A, (4-21) 
the dependence of V(kR , k c Ro,A) on /c c -Ro an d A is negligible, that is, if 

( k\ 2n+3 

P (1 \n;k/k c ^0,A^oo)^a n {-\ , (4-22) 

with a n some constant which depends only on the spectral index n. From the results of 
Section |3.1| , we see that this condition is only satisfied when n = —2. In this case we have 
(see Eq. (|3^35|) ): 

1 ^ 

V(kRo; n = -2) * — (LR )- 2 + — (/ci2 ) _1 , (4-23) 

which takes the self-similar form. For n = —1, 1-loop diagrams yield logarithmic fee- 
dependent corrections to self-similar scaling (see Eq. fl3-28[ )), 



1 16 

V(kR ,k c R ) « —(kRo)- 1 + — — (fci?o) 

Z7T 11U207T 



3068 + 2135 In I — 



(4-24) 



The results for n = 0, 1 show a stronger power law breaking of self-similarity (see Eqs. ( |3-19| ) 
and (EM2D ): 



V(kR ,k c R ) 



(kR 



61 (kRo) 2n+z 



2itT(^) 3157r(n + l) r 2 (^) V^c 



(4-25) 



where rj = — (n + 1) is an exponent which measures the deviation from self-similar scaling, 
generally known as the anomalous dimension in the theory of critical phenomena ( poldenfekl 
|1992| , [Barenblatt 1979| ). A visual summary of these results is given in the next subsection, 
where we compare 1-loop NLCPT to the universal scaling hypothesis first proposed by 
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Hamilton et al. (1991). For the variance and average two-point correlation function, self- 
similarity breaking means that, in the expected scaling region given by Eq. ( |4-17| ), the 1-loop 
coefficients sW and (see Eq. ( |2-2(J| ) ) are scale-dependent; in Sec. |^, we found this to 
be the case when n > — 1. Thus, for n > —1, we find generally that self-similarity is broken 
in perturbation theory by the first non- linear (1-loop) corrections to the power spectrum. 

Our conclusions about self-similarity in NLCPT differ substantially from those of Jain 
& Bertschinger (1995). The primary difference stems from the fact that they take the small 
scale cutoff in the linear power spectrum, k c , to be time dependent in comoving coordinates. 
They choose k c to scale as the inverse correlation length, i.e., they fix k c Ro ~ 1 throughout 
the evolution. Their rationale for this is plausible: since NLCPT is expected to break down 
at scales below the correlation length R , within the framework of perturbation theory one 
should perhaps restrict contributions to the power spectrum to scales larger than R . Since 
Ro depends on time, the cutoff of the linear power spectrum must be chosen to scale as l/Ro 
in order to satisfy this requirement at all times. Using this condition on k c Ro, they correctly 
conclude that the power spectrum evolves self-similarly (to any order in NLCPT), provided 
there are no infrared divergences. At 1-loop, we can see this directly from Eq. ( |4-2C| ) : if 
k c R = 1 and if p^ is convergent in the limit A — > oo, then the power spectrum depends 
only on kRo and therefore obeys self-similarity. However, while Jain & Bertschinger did 
show that leading infrared divergences cancel as loop momenta become vanishingly small, 
they did not compute loop corrections to the power spectrum explicitly. In fact, as we saw 
above, for n = —1, —2, p^ 1 ' develops a singularity at k = k c , and the limit A — > oo does not 
exist. This necessarily introduces a dependence of the 1-loop power spectrum on A = k c /e. 
As a result, even if k c Ro is held fixed, self- similarity is broken, unless e is also assumed to 
scale with R as k c does; we do not know of any arguments (analogous to that above) that 
would require e to scale this way. 

We take a different viewpoint on this question. As Jain & Bertschinger note, it is 
not surprising that choosing k c Ro »s 1 throughout the evolution leads to self-similarity, 
since one is identifying the cutoff scale with the correlation length. In fact, this choice 
can be characterized as artificially introducing self-similarity into the problem via the initial 
conditions. Here, we are rather interested in the question of whether self-similarity arises 
from the dynamical evolution over some range of scales far from the cutoffs (e.g., that 
given by Eq. ( [4-2 ip ) when the initial conditions are not exactly scale-free. The basis of 
perturbation theory is the assumption that the statistical evolution of the large-scale fields 
can be determined perturbatively even when the density and velocity fields are highly non- 
linear on the smallest scales. Therefore, when k c Ro > 1, one should include perturbative 
contributions from scales below Ro and see whether the results are sensible and/or self- 
similar, rather than continuously remove this power transfer by hand on the grounds that 
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the theory is expected to break down. This assumption of course will break down at late 
enough times and/or small enough scales, but it can be tested, and its range of validity 
quantified, by checking its predictions against numerical simulations. A fortiori, our choice 
of fixed comoving cutoffs is close to the conditions in (scale-free) N-body simulations, where 
infrared and ultraviolet cutoffs are imposed by the box size and the grid scale, both of which 
are fixed in comoving coordinates. [This is not strictly true in the later stages of evolution for 
high-resolution codes (such as P 3 M and tree codes), where the small-scale cutoff is eventually 
set by the mean interparticle distance; when the small-scale clustering stabilizes, this yields 
k c oc a, which, while not constant, has the opposite time dependence from the scaling needed 
for self-similarity, k c oc a~ 2 ^ n+3 \] 

We expect that in models with relatively more small-scale power, i.e., larger n, the basic 
assumption of NLCPT will have a smaller range of validity: the power being fed from small 
to large scales soon invalidates the perturbative approach. In fact, the stronger breaking of 
self-similarity we have found in NLCPT at larger n points to this conclusion, because scale- 
free N-body simulations exhibit self-similar scaling to very good accuracy for n > —(1 — 3). 
To assess this issue quantitatively, we now turn to a comparison of NLCPT predictions with 
the "universal" scaling observed in scale-free numerical simulations. 

4.4. Comparison with Simulations: the Universal Scaling Hypothesis 

The universal scaling hypothesis ([Hamilton et al. 1991| ) asserts that the non-linear 
evolution of the average two-point correlation function can be obtained from its linear coun- 
terpart by a universal scaling relation. This relation has been derived empirically from the 
study of numerical simulations for scale-free and CDM power spectra, but its main features 
can be understood on physical grounds ( [Nityananda fc Padmanabhan 199"4] ). These results 
have been extended by Peacock & Dodds (1994) to a universal relation between linear and 
non-linear power spectra and by Jain et al. (1995) to include the dependence of the scaling 
relation on the spectral index n. In its current version, the universal scaling hypothesis for 
the dimensionless power spectrum, 



A(Jfe) = 4nk 3 P(k), 



(4-26) 



relates the linear (A^) to the non-linear evolved (A e ) spectrum via 





(4-27) 
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where the simulations are empirically fit by setting B(n) = (1 + n/3) 1 ' 3 and 

(l + 0.6z + z 2 -0.22?-l.bz 7 ' 2 + * i \ 1/2 , Ano , 

$(z) = z . (4-28) 

V ; I 1 + 0.00372 3 j v ; 



The ansatz in Eq. (|4-27 ) is equivalent to the hypothesis that the mean dimensionless pairwise 



velocity is a universal function of the average two-point correlation function ( [Hamilton et| 
|al. 1991] , [Nityananda fc Padmanabhan 1994]) . Note that the linear spectrum is mapped 



to the non-linear spectrum at smaller scale (k e > k), reflecting the fact that non-linearly 
evolving perturbations shrink in comoving coordinates. The specific relation between k and 
k e in ( f4-27| ) can be obtained from the evolution equation for the average correlation function 
(conservation of pairs), according to which the mean number of neighbors within distance 
r e of a particle is time independent if r e decreases in time according to r 3 = r 3 [l + £ e (?"e)]; 
here r is a fixed scale, the radius of a sphere which encloses the same number of neighbors 
in the linear regime. 

We now compare these expressions with the perturbative calculations of the previous 
sections. We can write the dimensionless power spectrum up to 1-loop corrections as (see 
Eq. 



A e (^ ,M? ,A) = + ] W " + l kR °] p (1) (n; k/k c , A). (4-29) 



In Fig. [15] we compare the 1-loop NLCPT relation between linear and non-linear dimension- 
less power spectra with the universal scaling hypothesis (USH) for spectral index n = —2. 
The different dotted curves correspond to the predictions of NLCPT for k c Ro = 10°, with 
a = 0.5 — 2; since 3 < k c Ro <C A = 10 3 , cutoff effects in the linear evolution of the cor- 
relation length are negligible [see Eq. (|4-17 ) and Fig. (0)]. Each curve extends down to 



scales such that (kRo) max = 10 a-1 / 10 . The agreement between NLCPT and the USH (solid 
curve) is impressive all the way down to scales where A e w 10, far beyond the domain of 
linear perturbation theory (see also Lokas et al. (1995) for a comparison between NLCPT 
and N-body results for n = —2). For smaller scales than this, the 1-loop power spectrum 
overestimates the non-linear power spectrum. Note, however, that the turnover at A e w 10 
in the USH curve marks the onset of the transition to the regime described by the stable 
clustering hypothesis, and therefore disagreement with NLCPT beyond this point is not 
surprising. The NLCPT curve for a = 2 shows a slight disagreement with self-similarity 
at the large-scale end, which corresponds to the fact that, at this stage of evolution, these 
scales are not negligible compared to the size of the system. Otherwise, the agreement with 
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self-similarity is excellent. For comparison, we also present the perturbative results to 1-loop 
for the Zel'dovich approximation (ZA) for k c Ro = 100 (dot-dashed curve). The agreement 
between ZA and the USH is not very good, although the results of the N-body simulations 
of Jain et al. (1995) display a scatter comparable to the difference. In fact, their earlier 
(higher redshift) outputs seem to fit the ZA predictions shown here, suggesting that some 
of their scatter for n = —2 may be due to transients from the ZA initial conditions of the 
simulations. 



In Figs. ^ - [18], we show the same comparison for spectral indices n = —1,0,1. As 
expected from the discussion in the previous sections and the fact that the USH is self- 
similar by construction, for n > — 1 there is obvious disagreement between NCLPT and 
USH. This discord becomes more evident as time evolution proceeds (a change in the scale 
factor by Aa corresponds to Aa = lo(" +3 ) Aa / 3 ), as the scales considered decrease, or as the 
spectral index increases. The dotted curves that drop after reaching a maximum end close 
to the scale where the 1-loop corrections dominate over the tree-level contribution, driving 
A e negative (this happens at kR Q < 1). A similar conclusion regarding the disagreement 
between NLCPT and the USH has been recently noted by Bharadwaj (1995), who examined 
the hypothesis that the mean dimensionless pairwise velocity is a universal function of the 
average correlation function. From the analysis of the BBGKY hierarchy in the weakly 
non-linear regime, he showed that this hypothesis is not obeyed at the 1-loop level for linear 
power spectra Pn(k) oc kexp(-k), consistent with our results for n = 1. In fact, we find 
that non-linear perturbation theory fares worse than linear theory for n > — 1. 

A comparison of the USH with NLCPT to 1-loop has also been carried out by Lokas 
et al. (1995), who concluded, based on numerical calculations of the s^ 1 ' coefficients for 
Gaussian smoothing, that the agreement was good. Their main focus, however, was to 
understand whether non-linear evolution decreases or increases the growth rate of structure 
with respect to that predicted by linear perturbation theory. In this regard, NLCPT qualita- 
tively agrees with the USH in the sense that both predict that the growth rates are increased 
over linear theory when n = —2 and decreased when n > (with n = — 1 showing marginal 
behavior). However, as we have seen, the quantitative agreement is poor for n > — 1: the 
scale-dependence of the s' 1 ' coefficients in NLCPT disagrees with that of the USH, which by 
self-similarity requires ~ const, in the scaling regime given by Eq. ( 4-21 ). 



In fact, one can extract the s^ 1 ' coefficients predicted by the USH from Eqs. ( [4-2 7 ) 
and ( 4-28|) by using a small-A e expansion, 



A (D 



0.3 



n + 3 



1.3 



n + 3 



Ai, (4-30) 
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[compare to Eq. ( 4-29 )1 ■ Using Gaussian smoothing, this gives 



(USH) 



2 T(n + 3) 

r 2 (¥) 



0.3 



n 



1.3 



n + 3 



(4-31) 



which yields Sq* ~ 0.58, —0.32, —3.56, —13.52 for n = —2, —1, 0, 1 respectively. These results 
should be taken with caution, however, because the range over which the next-to-leading 
order term in the fitting formula ( |4-28| ) dominates is quite narrow ( |Jain et al. 1995| ) and 
therefore sensitive to uncertainties in the fitting formulae. To quantify this uncertainty, 
one can, for example, calculate x^h f° r n — ~ 2 from the result in Eq. ( |4-30| ) , which gives 



x 



(i) 

TH 



-2 



0.50; on the other hand, a more direct calculation from the fitting formula 



for the average correlation function (analogous to Eq. (|4-28|) ) given in |Jain et al. 1995 
yields x^(—2) m 0.75, in remarkable agreement with the perturbative result, Eq. fl3-40|) . 
For n — — 1, 0, 1, the USH fitting formula for the average correlation function yields x^h ~ 
—0.04, —0.55, —0.98. We also note that these results and Eq. ( |4-31| ) predict a similar behavior 
regarding the change of sign of the 1-loop corrections as n increases, as pointed out by Lokas 
et al. (1995). 



4-4- 1- Self- similarity via Dimensional Regularization 

In the results presented so far, we imposed infrared and ultraviolet cutoffs on the initial 
power spectrum to obtain finite integrals. For n = —2, however, we have seen that the 
resulting 1-loop power spectrum in the scaling regime is independent of the cutoffs and 
thus self-similar; this suggests that in this case the loop corrections can be alternatively 
calculated without the need to employ cutoffs. In fact, if the spectral index is in the range 
—3 < n < — 1, one can show that 1-loop NLCPT preserves self-similarity by considering 

in the no-cutoff limit and using dimensional regularization to regularize the required 
integrations (see Appendix B). Dimensional regularization is much simpler than the cutoff 
approach, because it obviates the need for complicated constraints on the angular integration 
variable [see Eq. fl3-5| )1- It also makes possible the analytic calculation for non-integer values 
of n. 

In Fig. [H], we show z^, derived via dimensional regularization for —3 < n < — 1, 
and compare it with the predictions of the USH taken directly from the fitting formula 
for the average correlation function (labeled USH - £) and from the dimensionless power 
spectrum ansatz (USH - A). We see that the agreement is very good over most of this 
range. The NLCPT results diverge as n — > — 1, as expected from the logarithmic divergence 



-34- 



found as k c — > oo in the cutoff calculation. As n — > —3, the NLCPT result goes to the 
unsmoothed value = s (1) = 4063/2205 « 1.843 (§ coccimarro fc Fneman 1996[) ; this can 



be understood from the fact that large scale fluctuations become dominant and averaging 
over small scales has a negligible effect. (By contrast, the USH has a singularity at n —>■ —3 
due to the way that the n-dependence is parametrized through B(n) in the USH fitting 
formula; since the USH was extracted from simulations with n > —2, it should probably 
only be trusted in that range.) The linear part of the NLCPT curve is well described by an 
expansion about n = — 3: 

oo) = ^1 - H| (n + 3) « 1.843 - 1.168 (n + 3). (4-32) 

Similar computations can be carried out in NLCPT for Gaussian smoothing and for the 
smoothed coefficient; the results show essentially the same features as displayed in Fig. [19 



5. Summary and Conclusions 

We have calculated the one-loop (first non-linear) corrections to the power spectrum, 
average two-point correlation function, and variance of the density field, including smoothing 
effects for top-hat and Gaussian filters, for scale-free Gaussian initial conditions. For the 
power spectrum, these results should replace the expressions given in Makino et al. (1992) 
for the P22 contribution on lengthscales smaller than the inverse ultraviolet cutoff of the linear 
spectrum (k c < k < 2k c ). Our results for the one-loop corrections to the variance extend 
those of Lokas et al. (1995) to top-hat smoothing and correctly go over to the unsmoothed 
values found in Scoccimarro & Frieman (1996) when the smoothing radius approaches zero. 

We found that, when formulated with "fixed" comoving cutoffs, non-linear perturbation 
theory beyond tree- level does not obey self-similar scaling when the spectral index n > — 1. 
As a consequence, in this spectral range, NLCPT disagrees with the results of N-body 
simulations embodied in the self-similar universal scaling ansatz (USH). For —3 < n < — 1, 
however, we found that one-loop corrections do obey self-similar scaling and are in excellent 
agreement with the USH down to lengthscales where the dimensionless power spectrum 
A e (k) ~ a 2 (R ~ 1/fc) ~ 10; below this scale, the N-body results make the transition to the 
highly non-perturbative regime described approximately by stable clustering. 

We interpret the breaking of self-similarity for n > — 1 as a signature of the breakdown 
of the fundamental assumption that the large-scale evolution of the density field can be cal- 
culated perturbatively when there are highly non-linear fluctuations on small scales. In this 
instance, self-similarity breaking can be thought of as arising from an ultraviolet divergence 
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in the 1-loop corrections as the small-scale cutoff k c of the linear power spectrum goes to 
infinity. This problem is more severe as the spectral index n increases, which is the expected 
behavior due to the increase in small-scale power. 

Since for n > — 1, the relative power on small scales is large, the question arises of 
whether the disagreement between the perturbative results and numerical simulations is 
due to the breakdown of perturbation theory (as we argued above) or to the use of the 
single-stream approximation, which neglects the effects of pressure gradients due to velocity 
dispersion. When these effects are included, one expects that the anisotropic velocity disper- 
sion (non-radial motions) associated with small-scale substructure can inhibit the collapse 
of large-scale perturbations — the "previrialization" effect ( pavis fc Peebles 1977| , [Evrard 5| 



Crone 1992| , [Lokas et al. 1995] , [Peebles 199(J| ). In fact, as n — > — 1 from below, s^ 1 ' changes 



sign and becomes negative, showing that such an effect is already present at some level in the 
single-stream approximation. In addition, the terms neglected in the single-stream approach 
have been included by Bharadwaj (1995), who nevertheless finds exactly the same results as 
in the single-stream approximation. Therefore, if velocity dispersion terms indeed become 
important, it appears that their effects cannot be treated perturbatively. This suggests that 
the interpretation given above is the correct one, i.e., for n > — 1 the non-linear evolution of 
the density power spectrum is inherently non-perturbative, and perturbative methods are of 
little use. 

Given the excellent agreement between 1-loop NLCPT and N-body results for n < — 1, 
one is naturally led to ask whether this agreement can be improved upon, i.e., extended 
to smaller lengthscales, by going to next order (2- loop corrections). We regard this to be 
unlikely, since the USH suggests that the physical picture involved in this small-scale regime 
is well described by stable clustering, which cannot be obtained perturbatively starting from 
linear solutions. Another interesting question is whether the validity of loop corrections in 
NLCPT for higher order cumulants (such as the bispectrum, skewness, etc.) will follow the 
same pattern with spectral index. We will address this issue in future work. 
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Eqns.([3^-|3^49|). 



-36- 



A. Expressions for I a and 1^ 

Recall the definitions given in the text (see Eqs. ( |3-37| )) : 



I a (n,z) = / u n+2 W 2 {u) du, 
Jo 



L(n, z) = / u n+2 W{u) du. 
Jo 



Here we tabulate these integrals as a function of filter and spectral index. 



A.l. Top-Hat Smoothing 



(Al) 
(A2) 



4(-M) 



Wz 5 

+4z 5 Si(2z) 
9 9 
4 + 8? 



3 - 5z 2 + (3 - z 2 + 2/) cos(2z) + (6z + z 3 ) sin(2z) 



1 - 2z 2 + cos(2^) + 2;zsin(2z) 



4z 2 ' 



+22sin(2^)], 



(A3a) 
(A3b) 



— [-1 - 3z 2 + (1 + z 2 ) cos(2;z) + 2zsin(2^) + 22 3 Si(2^)], (A3c) 
--(1 - 7 e - In 2) + -^[-1 + cos(2z) - 22 2 Si(2^) + 2z 2 hx{z) 



(A3d) 



where Si (z) and Gi(z) denote the sine and cosine integrals, defined by 



Si(z) 
Ci(z) 



sm{u) , 
du. 

11 u 

00 cos(u) 



du, 



u 



(A4) 
(A5) 



and 7 e ~ 0.577216... is the Euler-Mascheroni constant. As z — > 00, we have Ci(z) — > and 
Si(^) — > 7r/2. Using the following expression of a top-hat filter in terms of Bessel functions: 



W T u(u) = 3^| M - 3/2 J ; 



3/2 W, 



(A6) 
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we have (for any n in the range —3 < n < 1): 

in _ 2 r(l - n)T[(n + 3)/2] 



1,(71,00) = 9vr2 r 



T 2 (l -n/2)r[(5-n)/2] : 
where T(x) denotes the gamma function. For Ig(ji,z) we obtain: 



| cos(z) sin(z) 

01 W H : ~ 



z z^ 



It(-l,z) = 3[l-sin(z)/z], 

I((Q,z) = 3[Si(z) -sin(z)], 

J € (M) = 3 [2 - 2cos(z) - zsin(z)]. 
Using Eq. (|A|), we have (-3 < n < 0): 

r , n x o /-on r[(n + 3)/2] 

^,oo) = 3v ^2 r(1 _ n/2) - 
A. 2. Gaussian Smoothing 



(A7) 



(A8a) 
(A8b) 
(A8c) 
(A8d) 



(A9) 



U-2,Z) 


2 v ' 


(AlOa) 




= \il-eM-z 2 )} 


(AlOb) 




= ^ erf (2) - ^zexp(-z 2 ) 


(AlOc) 


Ia(hz) 




(AlOd) 
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where erf (z) denotes the error function, 

2 



erf(z) = — = I exp(-u 2 )du, (All) 

V7T J 

and for I^(n, z) we have: 

J f (n, ^) = 2 (n+3)/2 J CT (n, z/v^). (A12) 



When z — > oo we obtain: 



J ff (n l0 o) = ir(^V (A13) 



B. Dimensional Regularization 

To obtain the low- A; behavior of the power spectrum for n < — 1, one can use dimensional 
regularization to simplify considerably the calculations, and get one-loop coefficients such as 



x 



n; oo) for n in the range —3 < n < —1. Since we are interested in the limit k c — > oo, all 
the integrals run from to oo, and divergences are regulated by changing the dimensionality 
d of space: we take d = 3 + e and expand in e < 1. For power spectrum calculations, we 
need the following standard formula for dimensional-regularized integrals QUollins 198-4) ) : 



(g 2 )-[(k- q ; 



2]u 2 



T(d/2 - u 1 )T(d/2 - z/ 2 )r(z/ 1 + v 2 - d/2) 

r(vi)r(i/ 2 )r(d -v x - u 2 ) 



7T d/2 k d ' 2ui ~ 2 " 2 . (Bl) 



When using this equation, divergences appear as poles in the gamma functions, which can 
be handled by the following expansion (n = 0, 1, 2, . . . and e — > 0): 



r(-n + e) 



-l) n 



^ + ^(n + 1) + € - ( y + i?(n + 1) - i)\n + 1) J + 0(e 2 ) 



:b2) 



where r/>(:r) = d\nT(x)/dx and 



^(n + 1) 
^'(n + 1) 



1 1 
1 + - + ... + - 

2 n 

6 fr^ 2 ' 



7e 



(B3a) 
(B3b) 
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with ip(l) = — 7 e and if)'(l) = n 2 /6. We can write the one-loop power spectrum as (see 
Eq. (El) ): 



P (1) (k,T;n) = A 2 a\r) j d\ q n (2|k - q| n [F 2 (s) (k - q, q)] 2 + 6k n F 3 (s) (k, q, -q)) . (B4) 

In the numerators of the integrands, we can use the relation 2 q • k = — (k — q) 2 + k 2 + q 2 
to rewrite Eq. ( |B4| ) exclusively in terms of integrals of the form ((Bl|). The resulting 1-loop 
power spectrum contributions for —3 < n < — 1 are 



= / r(5/2-n)r 2 [(n-l)/2] 3r(3/2 - n)T[(n - l)/2]r[(n + l)/2] 
221 ' ' J " V 2r 2 (2-n/2)I>- 1) r(l - n/2)r(2 - n/2)I» 

29r(l/2 - n)r 2 [(n + l)/2] lir(l/2 - n)r[(n - l)/2]r[(n + 3)/2] 



4r 2 [l - n/2]r(n + 1) 4r(2 - n/2)r(-n/2)r(n + 1) 

15r(-l/2 - n)r[(n - l)/2]r[(n + 5)/2] + 15r(-l/2 - n)r[(n + l)/2] 



x 



2r(-l - n/2)r(2 - n/2)r(n + 2) 2r(l - n/2)r(-7i/2) 

r[(n + 3)/2] _ 25r(-3/2 - ra)r[(n + l)/2]r[(n + 5)/2] 25r(-3/2 - n) 
T(n + 2) r(-l - n/2)r(l - n/2)r(ri + 3) 4r(-2 - n/2) 

l)/2]r[(n + 7)/2] 75r(-3/2 - n)T 2 [(n + 3)/2] \ 7T 3 / 2 A 2 q 4 (r) 
X r(2-n/2)r(n + 3) + 4r 2 [-n/2]r(n + 3) / 49 
xk 2n+3 , (B5) 



= / r[(n + l)/2]r[(l - n)/2] 19r[-(n + 3)/2]r[(n + 5)/2] 

131 ' ' 1 \ 84r(l - n/2)r(l + n/2) 84r(-l - n/2)V(3 + n/2) 

r[-(n + 5)/2]r[(n + 7)/2] 5T[-(n + l)/2]r[(n + 3)/2] 

12r[-2 - n/2]r(4 + n/2) 28r(2 + n/2)r(-n/2) 
r[(n-l)/2]r[(3-n)/2]' 



7T 



42r(2-n/2)r(n/2) 
Note that this implies that oc /c 2ra+3 , as required by self-similarity 



A 2 a\r) k ln+A . (B6) 
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Fig. 6. — Same as Fig. but for n = —2 and A = k c /e = 100. For reasons of clarity, we do 
not show the behavior for x < 2A _1 (k < 2e), where P22 turns over to the scaling given in 
Eq. ( |3-7D . Similarly, P13 = when x < A -1 . 
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Fig. 7. — One-loop corrections to the average two-point correlation function x^ 1 ' and variance 
sW for n = —2 as a function of k c R, for top-hat (TH) and Gaussian (G) smoothing. Solid 
curve shows s^hi dotted curve corresponds to Sq\ dashed curve shows Xq\ and the long- 
dashed curve x^h- The dot-dashed line shows the large-scale approximation given in Eqs. 
( |3-39|) and (|3-41|) . For k c R = we recover the unsmoothed results s^ 1 ' = x^ ~ 1.82 
( |5coccimarro fc Frieman 1996|) . 
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Fig. 8. — Same as Fig. |7] but for n = — 1. Dot-dashed curves correspond to the large-scale 
approximation to the 1-loop coefficients, Eqs. ( |3-43| ) and ( |3-44j ) . We do not show x^, which 
undergoes large oscillations. 
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Fig. 9. — Same as Fig. |7| but for n — 0. Dot-dashed curves correspond to the large-scale 
approximation to the 1-loop coefficients, Eqs. ( |3-46| ) and ( |3-47| ). 
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Fig. 10. — Same as Fig. |7] but for n — 1. Dot-dashed curves correspond to the large-scale 
approximation to 1-loop coefficients, Eqs. ( |3-49| ) and ( |3-5C| ). 
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Fig. 11. — AI a ^(— 2, k c Ro, A) as a function of k c R for top-hat and Gaussian smoothing (see 
Eqs. ( ^4- 1 1| ) and ( |4-13| ) for definitions). Dashed lines represent the self-similar case (constant 
AI) in which the limit A = k c /e — > oo has been taken (e — > 0, k c — > oo). Dot-dashed curves 
denote the choice A = 10 3 , solid curves correspond to A = 10 2 . The differences between 
these curves and the A — > oo lines measure the deviations of the correlation length from 
self-similar evolution. 
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Fig. 12. — Same as Fig. [XT] but for n = 



-1. 
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Fig. 13. — Same as 
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11 but for n — 0. 
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Fig. 14. — Same as Fig. [IT] but for n—1 
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Fig. 15. — The evolved dimensionless power spectrum A e (k e ) as a function of the linear 
dimensionless power spectrum Ae(k) for n = —2. The scales k e and k are related according 
to Eq. Q4-271) . The solid curve shows the universal scaling hypothesis (USH); dotted curves 
correspond to the predictions of non-linear cosmological perturbation theory (NLCPT) in- 
cluding 1-loop corrections. The different NLCPT curves correspond to different values of 
k c Ro in the scaling region where self-similarity is expected to hold. The solid straight line 
corresponds to linear perturbation theory, whereas the dot-dashed curve is the prediction of 
the 1-loop Zel'dovich approximation for k c R = 100. 




Fig. 17. — Same as Fig. 



but for n = 0. 
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Fig. 19. — One-loop correction to the average correlation function, x^li, as a function of 
spectral index n in NLCPT (dotted curve) and in the USH obtained from the correlation 
function ansatz (solid curve labeled USH - £) and from the dimensionless power spectrum 
ansatz (USH - A). 



